Vortex melting and decoupling transitions in YBa 2 Cu 4 08 single crystals 
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The vortex correlation along the c-axis in high quality single crystals of YBa2Cu40s has been inves- 
tigated as a function of temperature T in different magnetic fields, using the quasi-flux transformer 
configuration. A simultaneous sharp drop associated with the vortex lattice melting is observed in 
both the primary and secondary voltages(Vt op and Vf> t)- Just above the melting temperature, the 
vortices form three-dimensional line liquid with the correlation length along the c direction L c < t, 
the sample thickness. The temperature where a resistive peak in Rbot develops corresponds to the 
decoupling temperature T<j at which the vortices loose their correlation along the c-direction and 
they dissolve into the two dimensional pancake vortices. The H-T phase diagram for the YBa2Cu408 
single crystal is obtained. 
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I. INTRODUCTION 

The vortex dynamics in the mixed state of high T c 
superconductors (HTSCs) remains a topic subject to in- 
tense investigation because of its importance for both the 
fundamental research and future applications. It is now 
well established that the thermal fluctuation induces a 
second order transition from a high temperature vortex 
liquid phase to a low temperature vortex glass phase for 
the vortex matter in a superconductor with a strong dis- 
order. In a clean superconductor the second order tran- 
sition is replaced by a first order transition from a vortex 
liquid to an Abrikosov vortex latticed 

Flux melting in various high quality high T c super- 
conducting single crystals h|aS|-.been observed by differ- 
ent experimental techniques It demonstrates itself as 
a sharp resistive drop or a jump in the magnetization. 
An important issue concerning the melting transition is 
that if the vortices lose their coherence in the c direc- 
tion during the melting transition or they will keep the 
c-axis correlations intact and then lose them at a still 
higher temperature a This issue has been widely pursued 
by doing transport measurements with the dc flux trans- 
former configuration— So far the results remain contro- 
versial. Doyle et al.E3 and Fuchs et alO found a coinci- 
dence of the melting transition and a decoupling one on 
Bi 2 SifepaCu208 (Bi22f 2)-single crystals. However, Wan 
et al.lia and Keener et al.t3 did similar measurements on 
Bi22I2 single crystals. They concluded that the melt- 
ing transition took place in a two-stage fashion, i.e., the 
rigid lattice at low temperature first melts into a three di- 



mensional (3D) liquid followed by a decoupling transition 
uppa an increase in the temperature. Recently Blasius et 
performed /x-spin rotation measurements on Bi2212 
single crystals with different oxygen contents and they 
obtained evidence for a two stage transition of the vortex 
matter as a function of temperature under equilibrium 
conditions. On the other hand, in the YBa2Cu307_a 
(YI23) system a coincidence of vortex melting and loss 
of vortex correlation along the c-axis was observed by 
Lopez et al.E-3 A decoupling transition was not observed 
in the optimally doped YI23 system since the vortex in 
this system is always in the 3D region. 

Underdoped YBa2Cu40s (Y124) is naturally stoichio- 
metric and untwinned, it has a moderate anisotropy pa- 
rameter 7 as compared to those of Y123 and Bi22I2. The 
weakness of the pinning strength in this material reduces 
the contribution from spurious effects in the melting tran- 
sition in the vortex system. Thus it is a good prototype 
material for the investigation of the interplay between 
the vortex interaction and dimensionality on the vortex 
dynamics. Recently, we observed a fkst order melting 
transition in the YI24 single crystals.E3 It is of great in- 
terest to extend our previous work to examine the vortex 
dynamics along the c-axis. Here we report the transport 
data obtained in simultaneous measurements for both the 
primary (Vtop) and secondary (Vbot) voltages using the 
dc transformer configuration, in order to investigate the 
vortex correlation along the c-axis in the YI24 system. 
It is demonstrated that in YI24 single crystals, vortex 
lattice (VL) melts into a 3D vortex liquid via a first or- 
der transitionand and that the decoupling transition is 
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seperated from the melting one. Our results support the 
two stage transition scenario that the vortices first un- 
dergo a melting transition, followed by a disappearance 
of correlation in the transverse direction at a higher tem- 
perature. 



II. EXPERIMENTAL AND RESULTS 

YBa2Cu40s single crystals were grown by the high- 
pressure flux method as described previouslyEH These 
single crystals were needle-like with typical dimensions 
of 1.2x0.4x0.01 mm 3 . The zero field critical transition 
temperature T c (0) of the single crystals was about 78 K. 
Two single crystals were used for these measurements. 
Each crystal was carefully cleaved to obtain optically flat 
surfaces with the c-axis normal to the sample surface. 
Gold wires were attached to the top and bottom surfaces 
of crystal by using Platinum epoxy. Then the crystal was 
heated in air at 100 °C for 1 hour, the resultant contact 
resistance was typically below 0.5 fi. The electrical con- 
tact geometry for the measurement is shown in the inset 
of Fig. 1. The secondary voltage contact pair (5 and 6) 
on the bottom ab face was placed directly beneath the 
top face primary voltage pair (2 and 3). The distance 
between each voltage pair was about 0.5 mm. The resis- 
tance were measured using a low frequency (17 Hz) ac 
lock- in technique with an excitation current of 0.1 mA 
injected from contact 1 to 4 while measuring simultane- 
ously the voltage drops V 2 3 and V 56 with the help of a 
Kcithlcy 228 scanner. The magnetic field was generated 
by a 15 Tesla Oxford superconducting magnet and was 
applied parallel to the c-axis of the crystal throughout 
the measurements. 
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FIG. 1. The logarithmic plot of the temperature depen- 
dence of R a b and R c at different applied magnetic fields. From 
right to left: 11=0.2, 0.5, 1, 2, 3,4, 6,7, 8, 10, 12 T. Inset: the 
electrical contact configuration for the quasi flux transformer 
measurements. 

Shown in Fig. 1 are the representative Rt op (T) and 

Fig. 1 X.G. Qiuetal. 



R-fcot(T) (defined as V23/I14 and Vse/lu, respectively) 
curves taken at different applied magnetic fields up to 
12 Telsa for one of the Y124 single crystals. Similar re- 
sults were obtained for the other crystal. The resistivity 
was measured by first cooling down the sample below T c 
and then collecting the data while warming up the sam- 
ple. While the applied magnetic fields induce a negligible 
depression in the superconducting onset temperature for 
the in plane superconductivity, they appear to consider- 
ably suppress the temperature where Rh ot deviates from 
its normal state value. A sharp jump in H top with a 
magnitude of R/R„ ~ 10% (R„ being the normal state 
resistance) is clearly observed at each magnetic field. We 
attribute it to the occurrence of the vortex melting tran- 
sition. We define the melting transition temperature T m 
as the one where a sharp peak in dR/dT exists. The tem- 
perature T m defined by Rt op (T) and Rb t(T) are equal, 
i.e., T^ P =T^*. This is demonstrated more clearly by 
plotting the derivatives of Rf 0p and Rb ot with respect to 
T together, as shown in Fig. 2a. Below the melting 
temperature, R top =Rf, t. Immediately above the melt- 
ing temperature T m , Rt op becomes larger than Rbot- As 
the temperature increases, both Rt op and R;, ot increase 
and the difference between R top and Utot becomes larger 
and larger. Finally near a characteristic temperature T^, 
Rfcot reaches its maximum and at the same temperature, 
a small kink develops in Rt op , as can be seen in a sep- 
arated plot in Fig. 2b. It is interesting to notice that 
above T TO , dRt op /dT decreases with the increasing tem- 
perature, reaches a minimum at T^ before a second local 
maximum appears at T p . At exactly the same tempera- 
ture T p , dRbot/dT shows a local minimal. This behavior 
persists for the dR/dT curves at all applied magnetic 
fields. By plotting T m and T^ as a function of the mag- 
netic field, we obtain a phase diagram for this Y124 sin- 
gle crystal as shown in Fig. 3. Previously we have found 
that T m (B) could be well described by the anisotropic 
Ginzburg-Landau theory and the melting line could be 
fitted by an empirical formula B m (T)=31.4(l-T/T c ) 1 ' 44 
with T c =78.6 KH The fitting result is shown as the 
solid line in Fig. 3. In the mean time, angular dependence 
of the melting transition gives for Lindermann cdterium 
cl=0.14 and the anisotropy parameter 7=12. 4.Ej 



III. DISCUSSIONS 

The dc flux transformer configuration is useful means 
for probing the dimensionality and longitudinal corre- 
lation of vortices in the mixed state. It measures the 
velocity correlation for vortices running along the CuO 
planes of Y124 since the voltage drops V23 and V56 (and 
thus R( 0p and R{, t) in the mixed state is induced by 
the motion of vortices whose velocity is determined by 
the integrated Lorentz force on the vortex segments of 
length l c piercing the top and bottom surfaces, respec- 
tively. The relation between Rt op and Rf, ot is determined 
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FIG. 2. (a) dR top (T)/dT vs. T and dR^ydT^s. T at 
an applied magnetic field of 1 T; (b) Rt op and Rf, ot at 1 T. 
The temperature at the cross point of the two solid lines is 
the crossover temperature between two TAFF regions. 

by the velocity of the vortices in the top and bottom 
surfaces. When the vortex correlation length along the 
c-axis is longer than the sample thickness, the voltage 
drops at the top and bottom surfaces would be of the 
same magnitude. The fact that immediately above the 
melting temperature R top is larger than TLbot suggests 
that the l c is smaller than the sample thickness t above 
the melting transition. 

Above T c in the normal state, R top is larger than 
Rbot because the sample geometrical aspect, electrode ar- 
rangement and anisotropy result in a nonuniform current 
distribution, i.e., more current will pass through the top 
surface than on the bottom surface. The fact that in 
the normal state the ratio Rtop/R*ot is nearly constant 
suggests that the anisotropy of Y124 is nearly tempera- 
ture independent in this temperature region. Upon en- 
tering the superconducting state, the superconductivity 
will make an easy path for electrons to run across the 
CuO planes and therefore reduces the anisotropy. As we 
will discuss below, the vortices at this temperature region 
is the 2D pancakes, the voltage drop will be mainly de- 
termined by the anisotropy. Since we are measuring the 
resistance by passing a constant current, the decrease in 
anisotropy means that more current will pass through 
the CuO layers and that the current flowing in the bot- 
tom layers will increase. The decreasing anisotropy with 
decreasing temperature is reflecting in the reduction of 
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FIG. 3. H-T phase diagram extracted from the trans- 
port measurements. Solid circles are the points for the 
melting line, open circles are those for the decoupling line. 
Solid line is the fitting result with the empirical relation 
B m (T)=31.4(l-T/78.6) 144 . Dashed line is the decoupling line 
fitted with B d (T)=1405.4(l-T/78.6)/T. t is sample thickness 
and s the distance between the neighboring CuO planes. 

Fig.3 X.G. Qiuetal. 

dRh t/dT as shown in Fig. 2a. 

In the present measurement configuration, K top is 
equivalent to the in-plane resistance R a {, commonly mea- 
sured in a four terminal measurement. Our previous work 
on the in-plane resistivity demonstrated that for R top (T), 
above the melting temperature in the liquid state, there 
are two distinct parts where lnR top shows linear depen- 
dence upon 1/T with different slopes, which is a typical 
characteristic for thermally activated flux flow (TAFF) 
behavior £j We have identified the temperature T<j where 
the crossover between the two TAFF regions occurs as 
the decoupling transition through an activation enecgy 
analysis. The decoupling temperature is described byu 

Bcr = 16^k B Tse\l bl 2 ' (1) 

where s is the interlayer spacing between the neighboring 
CuO planes. With a 7 value of 12.4 a satisfied fit is 
obtained by Eq. (1) for the decoupling line shown as the 
dashed line in Fig. 3. 

We notice that the temperature where Rf, t shows a 
maximum corresponds to exactly the same temperature 
where a crossover between two different TAFF regions oc- 
curs as we have identified before. Such a peak resembles 
that of the out-of-plane resistance R c often appeacs-|iH. 
single crystals with large anisotropy such as Bi2212.E3t2l 
However it can be seldom seen in single crystal with mod- 
erate anisotropy such as fully oxidized Y123. Although 
there is no consensus on the occurrence of the peak in 
Rbot and R c yet, one possible reason for it could be the 
tunneling of vortices over the c-axis barrier in the liquid 
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state following by the complete coupling of CuO planes. 
Recently, there have been considerable efforts devoted to 
the study of, mterlayer charge dynamics by far infrared 
reflectivitylErEl Studies on the interlayer charge dynam- 
ics with the electric field component E polarized along 
the c-axis reveals that above a characteristic tempera- 
ture T p i, the superconductor behaves like a poor metal 
or an ionic insulator. The carriers are confined inside 
the CuO planes and the c-axis resistivity is controlled 
by the interlayer tunneling. The c-axis coherence is sig- 
naled by the appearance of a reflectivity edge associated 
with the lasephson plasma from the intrinsic Josephson 
j unctions B3'Ej Above T p i, the electron transport along 
the c-axis is incoherent. Below T p i, coherence is built 
up for the interlayer carrier transport. Therefore, given 
the high transition temperature and weak interlayer mag- 
netic coupling between the pancake vortices in the neigh- 
boring CuO planes, T^ would quite possibly correspond 
to T p i (the Josephson plasma frequency) at which the 
phase coherence and thus the superconductivity is built 
up for the entire system along the c-axis. The suppres- 
sion of the onset temperature for c-axis superconductivity 
by an applied magnetic field is due to the breakdown of 
the phase coherence along the c-axis by a magnetic field. 
This ia consistent with the model proposedJay Briceno 
et al.Ej and later developed by Suzuki et al.tJ who sug- 
gested that phase fluctuation induced dissipation con- 
trolled the dissipation and could be well described by 
the Ambegaokar-Baratoff theory. 

Although the picture of vortex melting has been widely 
accepted, the mechanism of the vortex melting has yet 
not well understood. Recently, Nonomura et al.Ej car- 
ried out Monte Carlo simulations of the three dimen- 
sional frustrated XY model, in which the melting temper- 
ature was determined by the helicity modulus along the 
c-axis. They found that the melting transition is propa- 
gated through vortex entanglement. At the melting tem- 
perature, the percentage of entangled vortices abruptly 
changes. Upon increasing the temperature, the entangle- 
ment length becomes smaller than the sample thickness 
and decreases rapidly. The dissipation is governed by the 
cutting of vortices above the melting temperature. As a 
consequence of the entanglement mechanism of the FLL 
melting, T m is scaled by the inverse of the system size 
along the c axis. Our results can be explained by their 
simulation results very well since we indeed observed a 
melting transition into 3D line liquid with its correlation 
length l c smaller than the sample thickness t. 

The_s£ortex melting scenario has been questioned by 
MooreE3 who argued that the melting transition could 
just be a second order crossover of the vortex matter 
from a three dimensional behavior to a two dimensional 
one when the l c in the vortex liquid becomes compatible 
to the sample thickness t. Because l c grows very rapidly 
as the temperature is lowered, the crossover region ap- 
pears narrow enough to be misinterpreted as a first order 
melting transition. For a cutting length of the sample 
thickness t, the cutting temperature can be calculated 
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by substituting L with £«0.01 mm. Here (/>o=2.07x 10~ 7 
Gauss cm" 2 
is the flux quantum, Eq = &o m K /167r 2 A 2 fc (T), n being 
the Ginzburg-Landau parameter and A a &(0) « 2000A, 
the penetration depth in the ab plane. The obtained 
vortex cutting line lies slightly below the melting line 
B m (T). Although our result can also be qualitatively ex- 
plained with this argument, it is not clear why such a 
sharp crossover is not observed in thin films (strong pin- 
ning disorder) since in any sample with a finite thick- 
ness larger than the distance between neighboring CuO 
planes, a crossover temperature is expected. One reason 
could be that the relatively smaller numbers of pinning 
centers in the single crystals (weak pinning) than in the 
thin films result in a larger Larkin correlation length L p 
inside which an ordered VL can persist .Ej In the case 
of thin films, the transverse correlation length is smaller 
than the average lattice spacing, only a second order glass 
transition can be observed.CJ This in turn suggests that 
a regular VL with a relatively large radius is a necessity 
for the occurrence of a melting transition. This is con- 
sistent with the fact that artificially introduced pinning 
centers will induce a small Larkin length and smear out 
the melting transition. Further work on the study of the 
relatiohship between pinning strength and Larkin length 
need to be done in order to clarify this point. 

Our results suggest the following physical scenario: in 
a clean Y124 single crystal, below the melting tempera- 
ture T m , the vortices form a rigid VL with l c 3> t. As 
the superconductor is warmed up, at T m , VL melting 
occurs. At this stage, the vortices lose their translational 
coherence while keeping the correlations along the c-axis. 
Thus the vortices in the temperature region just above 
the melting transition are in a 3D line liquid state. At 
this temperature region, the vortex form a kind of entan- 
gled vortices with their entanglement length l c smaller 
than the sample thickness t but much larger than the 
distance s between the CuO planes. The dissipation is 
governed by the vortex cutting and recombination. As 
the temperature increases, the correlation length in the 
c-axis becomes smaller and smaller and the difference 
between R top and Hbot gets larger and larger. Finally 
at the temperature T^, inter-plane decoupling occurs (l c 
<s) and the vortices lose the coherence in the c direction 
and the vortices dissolve into 2D pancake vortices that 
move independently in the individual CuO planes. 



IV. CONCLUSIONS 

In conclusion, we have observed the melting transi- 
tion in Y124 single crystals. The vortices following the 
melting transition are in the 3D line liquid state with a 
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correlation length in the c- direction l c smaller than the 
sample thickness but larger than the distance between 
the CuO planes of Y124. The temperature at which a 
resistive peak exists at Rf, ot is found to correspond to the 
interlayer decoupling transition temperature. Above the 
decoupling temperature, the vortices lose their coherence 
in the c-dircction and the dissipation would be governed 
by tunneling of 2D vortices across the CuO planes. Our 
results also suggest that the interlayer decoupling transi- 
tion is a continuous crossover rather than a sharp tran- 
sition. 
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